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I. Introduction 

For the aim of analysis and synthesis in state space, an option for direct measurement or reconstruction of the state variables 
is needed. If such information is hard to be obtained even in case of availability of methods for direct measuring, it is often 
considered the use of information needed received through other channels. For the purpose of collecting the needed 
information for the state variables of the system a wide use of observers of full and reduced order, steady -state error free or 
with steady-state error, is natural. Beside the number of methods and approaches every practical realization is still a matter of 
interest, especially when if immeasurable disturbing influence is to be taken into consideration for analysis and synthesis. 


II. Formulation of the problem 

The object considered for analysis is continuous linear time-invariant system with structure presented in Fig.l [1,3], 



Figure 1. Structure of the plant under consideration 


The continuous linear time-invariant system is considered in particular dual-mass DC electromechanical plant. For the 
purpose of the synthesis of state observer the plant will be described and analyzed in the state space. Considering eventual 
technical realization there is a set of state variables chosen, which elements corresponds to real physical variables as follows: 
x 1 - armature current t a ; x 2 - rotor speed aq ; x 3 - engine elastic torque M 12 and x 4 - working/actuating machine speed a> 2 
(plant output). As a result and with accordance with the structure of the plant, presented above the following system of 
equations describes the state variables of the plant: 


*1 =^z-(k vc u-c<t>x 2 )-^x 1 

ia K a 1 a 

x 2 = f (cOxj - x 3 — M ;i ) 

J 1 

*3 = Cl2<>2 - X 4 ) 

x 4 = -0 3 -m, 2 ) 

J 2 


( 1 ) 


In regards to the basis considered the plant can be described in the state space with the following equations: 

x — Ax + bu + v 
y — c T x + du 


( 2 ) 
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where x is an //-dimensional state vector, u is an /--dimensional input vector, and y is an ///-dimensional output vector. A, b, c, 
d are the system matrices of appropriate sizes and v is the disturbance input matrix. 

Therefore matrices A, b , c, d and v of equation (2) are: 
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( 3 ) 


Once the linear time -invariant system, presented by equations (2), is transferred to discrete model the below equations 
describe the discrete system in the state space: 

x(k + 1) = A*x(k ) + b*u(k ) + v\k) 

y(k) = c* T x(k ) + d*u(k) ^ 


The following conditions should be noted for the plant under analysis: 

(i) Real measurable state variables to be considered are the armature current i a and the rotor speed a> 1 ; 

(ii) The load torque M l2 in the common case is considered as constant or unknown variable with small fluctuations; 

(iii) State variables M 12 (engine elastic torque ) and a > 2 (working/actuating machine speed ) are to be considered 
immeasurable. 

Even if a certain devices are being used for measurement of uq, it is preferable such devices not to be implemented. For 
observation and monitoring purposes is usually natural the implementation of a current measuring device (for the state 
variable t a ), thus in this case there is no need of engine speed measurement. 

One of the premises for granting correct and error free estimations of state variables is to ensure that all of the variables 
influencing the system are measured. In this case in particular appears an element, for which direct measurement is not 
possible - the load torque M l2 . There is a need that this unknown and immeasurable element is measured or recognized - one 
of the options is by identifying it (observe it) as a new state variable. In this regard the entire research is aiming synthesis of 
discrete state observer by one measurable state variable and identification of immeasurable disturbing influence. 

III. Solving the problem 

The first step in finding a solution of the problem defined is check performance for the state variables observability, which is 
only available in case of measurable armature current i a — x 1 . If c* T — [1 0 0 0] and y = the observability matrix 

of the plant should be one of full rank: rank(Q* ob ) = n, where n is the plant order. 

As it was clarified above there is occurrence of disturbing influence for which direct information is not available and it is also 
considered immeasurable - load torque M l2 . Under those considerations the estimation of the state variables through the so 
defined observer will only be possible through feedback based on integral of the difference between the measurable state 
variables (in this case the measurable state variable is the armature current only x 1 — i a ) and their estimated values (x 1 ) [4], 
The feedback defined as stated above will be implemented for reconstruction of M l2 — * M l2 (fc + 1) = M l2 (k) + 
kobTqiyobQi) c llx(k)), where is coefficient of the feedback based on integral of the difference between the 
measurable state variable x 1 and its estimated value x 1 , and T q is discretization period, used for the discrete model of the 
linear time-invariant system under investigation. 

The observer equations then might be written as follows: 

x(k + 1) = A*x(k) + B* ob u ob (k ) + J*c*J b {x(k ) - x(k )) + v ob (k ) 

y ob (fe) = ell x(k) 


(5) 
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where c*J b = [l 0 0 0 ]; v ob (fc) = /*M, 2 (fc); u o6 (fc) = 


u(k ) 
Mn(fc) 


and f* denotes feedback of the state variables; f* is 


feedback based on integral of the difference between the measurable state variable x 1 and its estimated value x 1 ; M l2 (k) 
estimation of the disturbance influence M ;2 (k). 


If it is considered M i2 (fc + 1) = M, 2 (k) — k^TqX-^Qc) , where x 1 (k) is an error of the estimation of state variable Xj_, then 
the synthesized observer can be presented with the following equivalent matrix A* obequiv — A* —f*c* 0 T b . Defined in this 
manner the synthesis of the observer will be formulated as synthesis of the feedback matrix f* elements and k* ob . 


There is a requirement for the discrete observer to have zeroes poles, which will grant the shortest possible time for 
reconstruction of the state variables and identification of M ;2 (fc). The procedure for pole placement in predefined area is 
explained and simulated in details in [l].The same procedure is applied for synthesis of the above described observer. 

In this case the procedure is significantly simplified as the desired polynomial is the following: h* d {z ) = z 5 , e.g. the 
coefficients before grades of z which are less than 5, are zeroes. The elements of f* and k* ob could be obtained, considering 
that all the coefficients before the polynomial elements of the polynomial h(z ) = [zE — A * obequiv ] (where E is an identity 
matrix) are required to be zeroes, with exception the one which is before the highest grade, which is equal to 1. 


IV. Experimental results 


Dual-mass DC electromechanical plant shown in fig.l is operative through DC motor with the following parameters [1,3]: 


U nom — 220T; P nom — 0,3 kW\ n nom — 1000 tr/mirv, I anom = 2 A;u -» 0 -r 10F; 

R a = 20,8828ft; T a = 0,0126s; cO = l,702Fs; c 12 = l,702Abn; 

k pc = 22; J 1 = 0,042 kgm 2 ; ] 2 = 0,021 kgm 2 ; M (1 = 0,3404ATn; M l2 = 3,0636Abn 


For continuous -time state space system model we calculate: 
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Discrete-time state space model (4) could be obtained from the continuous one (2), considering (1) and (3), using MATLAB 
function for continuous to discrete model transition - c2d(.). Discretization period is defined T q = 0,1s which is chosen 
based on requirement T q < (0,05 -t- 0,l)t c . In this case t c (forecasted control time) is set to be equal to the expectation for 
the continuous-time closed-loop system [2], Considering these assumptions for the discrete-time state-space model (4) is 
calculated: 
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For the matrix of observability according to the chosen discretization period (T q — 0,1s) is calculated: 
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Q* ob 


1 -0,0291 

0 0,3313 

0 0,0569 

0 0,1383 


-0,0110 -0,0045' 

0,1129 0,0571 

0,0453 -0,0236 

0,4157 0,4751 . 


It is of full rank - rank(Q* ob ) — 4. Thus we can conclude that the discrete object is fully observable by one measurable state 
variable. 


For the observer equations (5) the following expression will be valid: 
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c* 0 T b = [ 1 0 0 0],v oh (k) = /*M/ 2 (fc) = [0,01675 
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In this case in particular if we have to meet the zero poles observer requirement the procedure for pole placement [1] within 
elliptical area will produce the following: 


- elongation of the center of the ellipse from the origin -> /r = 0; 


- major semi -axis -> = 0; 


- minor semi-axis -> r 2 = 0. 

With the defined elliptical area is ensured reconstruction of state variables and identification of disturbance influence M ;2 (/c) 
for t = nT q — 0,5s. 

As result of described synthesis procedure are calculated the following: 

k* ob = 79,7643; /* = [2,6558 -35,0648 12,416 - 29,928S] T 

In fig. 2 and Fig. 3 are shown immeasurable disturbance M i2 = 3,0636/Vm (applied at the 8 th second of the experiment) and 
the start processes of the state variables under signal value u = 9V, as well as their estimations. All the estimated values in 
the MATLAB environment and the plotted graphics are indexed - “est”. 



Figure 2. Immeasurable disturbance M 12 and its estimated value M 12 
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Figure 3. State variables and their estimations during the start processes 

In fig. 4 is presented the block-diagram of the system “plant-observer”, used as a basis for the simulating model built in 
MATLAB/Simulink for the purpose of experiment. 



Figure 4. Block diagram of the system “plant-observer” 
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V. Conclusions 

Analysis of the start processes shows that in 0,5s time are reconstructed the state space variables and identification of 
disturbance is performed. The steady state error is zero - practically no steady state error occurs. The synthesized discrete 
observer reconstructs the state variables and identifies the disturbance M l2 in terms of shortest period of time possible. 

Future developments 

There are number of methods applicable for reconstruction of immeasurable disturbing influence and further research will be 
focused on factorial design with further implementation in modal control algorithm. 
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